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Abstract 


A computational plasma aerodynamics model is developed to study the performance of an 
experimental laser propelled lightcraft. The computational methodology is based on a time- 
accurate, three-dimensional, finite-difference, chemically reacting, unstructured grid, pressure- 
based formulation. The underlying physics are added and tested systematically using a building- 
block approach. The physics modeled include non-equilibrium thermodynamics, non- 
equilibrium air-plasma finite-rate kinetics, specular ray tracing, laser beam energy absorption and 
refraction by plasma, non-equilibrium plasma radiation, and plasma resonance. A series of 
transient computations are performed at several laser pulse energy levels and the simulated 
physics are discussed and compared with those of tests and literatures. The predicted coupling 
coefficients for the lightcraft compared reasonably well with those of tests conducted on a 
pendulum apparatus. 
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Nomenclature 


Aj reaction rate constant of the j-th reaction 
Ap,Am matrix coefficients of transport equations 
c speed of light in a vacuum. 

D species diffusivity 

€v vibrational energy 

equilibrium vibrational energy 
F flux vector 

g, Gibbs’ free energy of species i 

H total enthalpy 

I intensity 

J emission coefficient 

kb Boltzmann’s constant 

kc Coulomb constant 

kv thermal conductivity of the vibrational energy equation 

Mw molecular weight 

me electron mass 

Ni quantum level population 

n refraction index 

7i unit normal vector 

He electron number density 

Pr Prantdl number 
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p gas static pressure 

pe partial electron temperature 

Qg,- elastic collision energy transfer 

Qr net radiative heat flux 

qe electron charge 

R gas constant 

7 position vector 

S source terms 

T gas temperature, K 

Te electron temperature, K 

Ti translational temperature, K 

Tv vibrational temperature, K 

Ui velocity components 

V velocity magnitude 

X, Cartesian coordinates 

chemical species mass fraction 
Greek Symbols 
At time step size 

eo permittivity 

(p dependent variable 

r control surface 


Kl 


laser absorption coefficient. 



A- 

Ai 

V,j 

Q. 

cb. 

COp 

p 

or 

nr 

Tij 


electron species thermal conductivity 

laser wavelength 

viscosity 

stoichiomatric coefficient 
control volume domain 
chemical reaction source term 
plasma frequency 
density 

Schmidt number 

vibrational-translational relaxation time scale 
shear stress tensor 


Introduction 

Beamed energy propulsion was first promoted by Kantrowitz . Since then, a propulsion 
system supported by a laser-sustained plasma has been the subject of many studies. i ne 

main advantage gained by laser propulsion is the low-weight system from decoupling the energy 
source from the vehicle, and high specific impulse from low fuel consumption. The first ground 
and flight tests of a ground-based laser propelled vehicle was reported in 1998,^’ in which spin- 
stabilized free flights with active tracking and beam control were accomplished to altitude of 4 
m. Six months later, the spin-stabilized vertical free flights reached 30 m."^^ These axisymmetric 
vehicles are of a special design among many innovative concepts in which a nosecone shaped 
forebody, an annular shroud, and a parabolic afterbody are the only major components. Vehicles 
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of such a design is hereby named as the “Laser Lightcraft” in this study. The layout of a 
computational grid for a Laser Lightcraft is shown in Fig. 1. The parabolic afterbody engine 
serves both as an aerospike nozzle and also the main receiving optic, whereas the annular shroud 
surrounds the ring focus of the parabolic optic. The Laser Lightcraft uses a Rocket-Based 
Combined Cycle (RBCC) to operate in both the rocket and air-breathing modes. When operating 
in an air-breathing mode, the specific impulse is infinity since no fuel is consumed. These 
successful tests'^' '^^ demonstrated the feasibility of launching small payload with a Laser 
Lightcraft. 

When focused beam energy strikes the shroud of a Laser Lightcraft, free electrons form to 
invoke Inverse Bremsstrahlung (IB) and the optical breakdown ensues. In order to achieve 
maximum intensity and to avoid plasma resonance, repetitive short pulses are used. As such, 
repetitive Pulsed Detonation Waves (PDW) are generated to propel the vehicle, hence the Laser 
Lightcraft engine is also know as a PDW engine. It can be seen that the underlying physics 
involved such as IB and PDW are so complicated that a simple system model is not capable of 
describing the phenomena, let alone predicting the thrust performance and the associated flow 
and thermal environments necessary for designing for the structure integrity. A detailed 
computational plasma aerodynamics model is therefore needed for the design and scaling of 
future Laser Lightcrafts. 

From May to July, 1999, performance data in the form of coupling coefficients were taken in a 
series of indoor tests for the Model #200 series vehicles. These test were performed under a 
cooperative effort among Rensselaer Polytechnic Institute, Air Force Research Laboratory and 
NASA Marshall Space Flight Center (MSFC), to enhance Laser Lightcraft propulsion research. 
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These data are the first of a kind benchmark quality data for computational model development. 
Armed with these data for comparison purpose, a computational plasma aerodynamics 
methodology is developed as a first attempt to predicting the thrust performance of a PDW 
engine propelled Laser Lightcraft. 

Computational Plasma Aerodynamics Development 

Some attempts have been made to analytically or computationally model the physics inside 
continuous wave (CW) devices. An early effort by Raizer assumed constant-pressure one- 
dimensional (ID) flow in air with heat addition by a laser. Thermal conduction was considered 
dominant and the radiation loss from plasma was ignored. Kemp and Root later extended this 
ID analysis to hydrogen and to include thermal radiation. Molvik, Choi, and Merkle extended 
this problem to a two-dimensional (2D) structured-grid hydrogen flow and implemented real ray 
tracing, but ignored thermal radiation. Jeng and Keefer^* did similar analysis and added thermal 
radiation. In this work, a multidimensional unstructured-grid computational plasma 
aerodynamics methodology is developed, with emphases on nonequilibrium effects and laser- 
induced physics such as laser - plasma interactions. Realistic laser absorption coefficients are 
used. The governing equations will be described first, with the auxiliary equations that 
computing the nonequilibrium and laser induced physics to follow. 

The Governing Equations 

To properly describe the plasma aerodynamics involved in laser propelled propulsion, the 
time-varying transport equations of continuity, species continuity, momentum, global energy 
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(total enthalpy) and electron energy are formulated and written in a Cartesian tensor form; 



where the shear stress vij can be expressed as: 




^ d Uj d Uj 



The beam energy absorption, beam energy transfer, and plasma radiation are invoked through the 
source terms of the global energy and the electron energy equations, whereas the plasma species 
are generated through the source terms of the species continuity equations. 
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Numerical Scheme 


The cell-centered scheme is employed that the volume surfaces are represented by the grid cell 
surfaces. The transport equations can be written in integral form as: 

F ndr=j^S^dQ. 

ot 


where q is the domain of interest and T denotes the surrounding surfaces; ^ is a unit normal 



vector of r in the outward direction. The flux function F contains the inviscid and the viscous 

flux vectors, And, ~ 

where V (f)^ is interpolated from the neighbor 
cells E and P. 

F=fN(t>-p^y(t> 


For the face e between control volumes P and E, the diffusive flux can be approximated as; 





f 

• n- 

V 





J 
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The finite volume formulation of flux integral can be evaluated by the summation of the flux 
vectors over each face. 



■ndT= 

J = k{i) 


where k(i) is a list of faces of cell i, Fij represents convection and diffusion fluxes through the 
interface between cell i and j, and ^Fj is the cell-face area. 

The convective flux is evaluated through the upwind-cell quantity by a linear reconstruction 

procedure to achieve second order accuracy: ~ ) where the subscript 

u represents the upwind cell and \f/e is a limiter used to ensure that the reconstruction does not 
introduce local extrema. The limiter proposed by Barth is used here. Defining 

^max ^ ‘Pmin = tnin((p,^ ,(/> j) (and assuming ipg is computed with ^^=1) the scalar 

y/e associated with the gradient at cell u due to edge e is: 


m 1 , 
m HI [ 1 , 



if 0>o 


if 0'\- <p <Q 
,/ 0 = 0 


Solution Procedures 

A general implicit discretized time-marching scheme for the transport equations can be written 
as below, 
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A/ 


* + 




NB 

= 1A,, 

nj=l 


vi+i ^ 


where NB means the neighbor cells of cell P. The high order differencing term and cross 
diffusion term are treated using known quantities and retained in the source term and updated 
explicitly. The second term of the right-hand side (RHS) is a perturbation which comes from the 
RHS. A predictor/corrector solution algorithm is employed to provide coupling of the governing 
equations. The discretized finite-volume equations form a set of linear algebra equations, which 
are non-symmetric matrix system with arbitrary sparsity patterns. The preconditioned Bi- 
CGSTAB''' and GMRES(m)^^ matrix solvers are used to efficiently solve the linear algebra 
equations. 


Auxiliary Equations 

Thermal Non-Equilibrium Energy Equations 

For high temperature flows, thermal non-equilibrium state may be important. In Landau and 
Teller’s derivation,^' a master equation is employed to describe the evolution of the population 
of quantum level A,. This master equation is written as: 

^ = w'z Kj^,Nj - n't K,^jN,-j = 0,1,2 /,„„ 

dt j=o j=o 

Results from the quantum mechanical solution of the harmonic oscillator are used to relate the 
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various quantum transition rates to one another, and then the master equation may be summed 
over all quantum states to obtain the Landau-Teller equation; 


Dt 


_d_ 

dxj 






An empirical expression (given below) is used to model the Landau-Teller relaxation time 
scale. First, the thermodynamic data for each species at every grid point are calculated using the 
following equations: 


R 

h 

RT 




= fl, + + 


T 


^ = a, \nT + aJ+^a,r- +^ajU^a,r +a, 


and the Landau-Teller expression for the vibrational relaxation time is calculated using the 6- 
coefficient specified. 


bj’'^ exp| 


[b.lT)' 


p(\-b^&\^-bJT)) 


where bi=7.1E-4, b2=0, b3=1.91E6, b4=0.3333, bs=0, b6=0 are used here. To solve this 
vibrational energy equation, the source term is linearized to result in an explicit term and an 
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implicit term. This treatment is important for an unconditionally stable solution of the equation. 
That is. 


Explicit source term = h 

^ LT 

P 

Implicit source term = “ 

^LT 

The vibrational temperature is used to influence the reaction rates of chemical reactions by 
assuming that the rate coefficients for dissociation are functions of the geometrical mean 

PI 

temperature between T and Tv. 

Non-Equilibrium Air Chemistry 

A general system of chemical reactions can be written in terms of its stoichiometric 
coefficients (vij and vij’ ) and the f' chemical species name (M/) of they"' reaction as 




The net rate of change in the molar concentration of species / due to reactions j, Xij, and the 


species production rate can be written as: 
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M,,' 


d)^ - Mw, ^ X^j 

j 


The forward (Arrhenius law) and backward reaction rates for each reaction is given by; 

/( E, 

Kf, -AT ‘ Exp( —) 

J) I DT 


Kb 


Ke, 


where Kej is the equilibrium coefficient 




Ke, = (RT) ‘ 


Exp 


I (V.. - V..) g. 

, = 1 ‘J 'J ‘ 

RT 


A point-implicit (operator splitting) method is employed to solve the chemistry system. For 
the breakdown of air. Park’s multitemperature air chemistry^* is based in this study. As shown in 
Appendix 1, this mechanism composes of the dissociation, NO exchange, associative ionization, 
charge exchange, electron-impact ionization, and radiative recombination reactions. With this 
mechanism, electrons are produced first by the associative ionization process and more electrons 
are produced by the electron-impact reactions. Since the number of electrons doubles in each 
such event, electron density increases exponentially, in the form of an avalanche. This 
mechanism provides the initial electron density for ignition. The combination of the avalanche 
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of electrons (produced from this mechanism) and the absorption of incident laser energy drives 
the subsequent optical breakdown. It can be seen that IB is embedded in the dissociation, 
associative ionization, and electron-impact ionization reactions. The reaction rates were 
validated*"' with experiments of postshock temperatures ranging from 20,000 - 60,000 k which 
are inline with the maximum computed breakdown temperatures. Notice heavy gas temperature 
T is used in place of Ta ( = for dissociation reactions in current study. 

Plasma Initiation 

Understanding the mechanisms responsible for plasma initiation (ignition) has been the subject 
of ongoing research. A spark gap, an extremely intense pulse of laser energy striking the focal 
point, seed molecules or particles, or a retractable tungsten target placed at the focal point “ have 
been used as the source of free electrons for initiation.*^' It has also been shown that plasma can 
be ignited quite easily off metal surfaces.*^’' Likewise in numerical modeling, initial free 
electrons are required for plasma ignition. This means determining an initial electron density or a 
threshold breakdown intensity. For example, Ergun*^^ described an initial electron density is 
necessary for modeling the ignition of a hydrogen plasma initiated with a retractable tungsten 
wire. In this study, a threshold breakdown intensity approach is devised in anticipation that 
different threshold breakdown intensity may be encountered for different surface material, target 
approach, and laser power. The advantage of this approach is that threshold breakdown intensity 
can either be calibrated or measured. For example, the effect of seeds on laser breakdown 
intensity was reported.*^' This is accomplished by performing parametric studies on plasma 
initiation with spark time, spark area, and spark power. A “spark” provides a fixed amount of 
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energy (spark energy) in an area (spark area) centered around the focal point for a fixed amount 
of time (spark time). When the strength of a spark reaches a threshold, or the threshold 
breakdown intensity is satisfied, a laser supported combustion (LSC) or detonation (LSD) is 
initiated. 

That means enough seed electrons are produced through the associative ionization reactions 
and the subsequent optical breakdown (ignition) is thereby possible. It is found that an ignition is 
sustainable when the energy absorbed by plasma reaches 15% of that of the delivered. Spark 
times computed based on that criterion are ranging from 0.4 ps to 1.2 ps. As expected, lower 
laser pulse energies require longer spark times for ignition. It is also reasonable to assume that 
the plasma initiates in an (spark) area enclosed by the two outer laser rays and an arc with an 
origin at the focal point. This is a cone (in a two-dimensional sense) and the (spark) radius is 
determined empirically as 1 mm. Finally, the spark power is also empirically determined as 15% 
of that of the incident laser - its theoretical maximum. 

Laser Radiation 

Geometric optics is used to simulate the local intensity of the laser beam, which is split into 
a number of individual rays. In the presence of absorption, the local intensity of each ray follows 
the Beer’s law: 


Electron is the only plasma species that absorbs the laser energy. The absorption coefficient of 
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the CO2 laser radiation, corrected for stimulated emission in the single ionization range, is 


approximated by the formula: 


R1 


5.72p^'ln[27(r/10‘*)^^^/7/''^] 

(y,/104)7/2 


A ray may change its propagating direction due to the inhomogeneous refractive index within 

El 

the hot plasma. The index of refraction is taken from Edwards and Fleck: 


1 - 


m C^TT 

- 



The refracted angle is associated with the refractive index through the Snell’s law:'^^ 
n, sin = «2 sin ^2 

where m and n2 represent the refractive indexes for two different control volumes, and 0 i and 02 
are the incident and refracted directions with respect to the normal direction of the interface 
between two volumes. 
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Non-equilibrium Plasma Radiation 

Treatment of radiative heat transfer of plasma is different from that of laser and the solution of 
the radiative transfer equation (RTE) is required. Neglecting transients and assuming a 
noscattering medium, the complete RTE becomes 


dI{s,Q) 

ds 


+ Kl{s,0.) ^ j(s) 


In equilibrium gases, the electronic energy-level populations are determined as a function 
of a uniquely defined equilibrium temperature according to a Boltzmann distribution, and k andy 
are related according to Kirchhoffs law as j(s) = fcl’’ (s) , where l\s) is blackbody intensity 

determined by Planck function. In non-equilibrium flow, however, all these simple relations no 
longer apply, whereas the non-equilibrium absorption and emission coefficients must be 
determined by the non-equilibrium populations of each energy level and transitions among 
various energy levels. 

Air plasma includes following atomic and molecular species: O, N, NO, N^, O^, NO^, N 2 , 
O 2 , N 2 ^, 02^, and e‘ which all contribute to the non-equilibrium radiative heat transfer. To 
determine k and j for these species, the following four radiative transitions must be taken into 
account: atomic line transitions, atomic bound-free transitions, atomic free-free transitions, and 
molecular transitions. Currently, there are two codes available which provide detailed 
information on these transitions. One is NEQAIR developed by Park’’^ and another is LORAN 
developed by Hartung.”' Because NEQAIR code involves large database, they are rarely applied 
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in practical problems. In this study, LORAN code is used to calculate the air plasma radiative 
properties. 

With the determination of k and j, the RTE can be solved by either deterministic or 
stochastic approach. Unlike the computational plasma aerodynamics governing equations, RTE 
is an integral differential equation and numerical treatment is, thus, different from those for 
differential equations. Currently, there are several methods available for solving non-equilibrium 
radiative heat transfer which include Monte Carlo method, P-1 method, quadromoment method, 
discrete ordinates method (DOM), discrete transfer method (DTM), etc. The Monte Carlo is 
accurate but too time-consuming for practical applications. The P-1 and quadromoment methods 
are only accurate for optically thick medium. The DOM and DTM are mathematically simple and 
can provide accurate results for all optical ranges if the discrete direction number is reasonable 
large. Thus, the DOM and DTM are very suitable for modeling radiation with a participating 
medium. For the laser propelled vehicle application, the high temperature region where non- 
equilibrium radiation is prominent is usually relatively small compared to the flow region 
considered. If the DOM is applied, RTE must be solved in the entire flow region, thus, much of 
the computational time is wasted in the regions where radiation is not important. However, with 
the use of the DTM, RTE can be solved in a designated region. Therefore, the efficient DTM is 
selected to solve the RTE. Detailed numerics of DTM is described in Ref. Cl. 

Plasma Resonance 

The plasma frequency is described as a property of a space-charge-neutral plasma by 
which the motion of the electrons in specific electrostatic oscillations is characterized. The 
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plasma frequency tOp of the electrostatic oscillation of the electrons is defined as®' 

2 2 , _ 

cop — HeCJe / nieEO 

It is seen that plasma frequency is a function only of the electron density. The electrostatic 
oscillations cause a resonance of the plasma to incident (laser) electromagnetic waves of the 
same frequency which are then totally reflected."^ A critical electron density above which the 
laser beam is totally reflected can then be obtained by equating the frequency of the incident laser 
with that of the plasma. 

Once the local electron density reaches the critical value, a “total reflection” condition has to 
be applied to the incident beam since the plasma becomes “opaque” and no beam can penetrate 
through the plasma front. This is achieved by assuming a total absorption at the plasma front and 
then a redistribution of the energy. Zhu^‘ noted that during the critical plasma resonance, the 
incident laser energy is to exciting the electron to a higher energy level, to reflecting and 
radiating diffusely, and to a lesser extent, escaping in the form of light by returning the excited 
electron to a lower energy level. The first energy redistribution mechanism can be viewed as an 
energy conversion, whereas the second and third ones can be grouped as an energy loss. 
Although the percentage of each is not theoretically known, an energy conversion percentage 
during the critical plasma resonance can be calibrated with the test results. In this study, a 30% 
energy conversion is used. 


Experimental Setup 
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The 10 kW Pulsed Laser Vulnerability Test System (PLVTS) CO 2 laser was used to provide 
the beam energy for Laser Lightcraft thrust experiments. The CO 2 laser delivered up to 400 J 
single-pulses at a pulse width of 18)as. Several variations of the basic Laser Lightcraft design 
(Model 200 series) similar to those described in Ref. M2 were examined and the test results of 
the bare aluminum Model #200-3/4 vehicle is chosen for this study. The impulse measurements 
were conducted with a pendulum apparatus. This technique employed a velocimeter coil, which 
was used in previous work to determine the impulse imparted to a flat plate using the PHAROS 
III laser at Naval Research Laboratory.^^ The Laser Lightcraft was suspended and suitably 
weighed before being subjected to a single pulse of energy from the PLVTS laser. The impulse 
hammer had been calibrated with clay balls and the data agree with each other well. Details of 
the laser, measurement technique, and pendulum apparatus can be found in Refs. Ml and M2. 


Computational Grid Generation 

Figure 1 shows the layout of a computational grid. As described in Refs. Ml and M2, the 
Laser Light vehicle consists of a forebody (nose), an annular shroud, and a PDW engine 
(parabolic optic). Only half of the grid shown in Fig. 1 is actually solved due to the 
axisymmetric formulation. An eight-zone structured grid was generated first using a grid 
generation code UMESH.'^^ The grid is then fed to an auxiliary program such that a one-zone 
unstructured-grid is resulted for the actual computation. Quadrilateral elements are used such 
that wall boundary layer can be computed if needed. High grid density is used in the inner 
shroud region for capturing the optical breakdown and pulsed detonation wave processes. 
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Results and Discussion 


A series of computations have been performed for average laser pulse energies of 75, 100, 
150, 200, 300 and 400 J. Maximum computed temperatures ranging from 51,000 to 60,000 K 
fall within the computed post-shock temperatures of 20,000 to 60,000 K using Park’s multi- 
temperature air chemistry.*^* Laser supported detonation occurs in all cases with maximum Mach 
number reaching 1.6. Laser supported combustion then ensues as the plasma front expands, 
whereas the Mach number decreases to subsonic level. Figure 2 shows the computed 
temperature contours and laser beam traces at an elapse time of 0.5 p,s and average laser pulse 
energy of 400J. It can be seen that the laser beam reflects specularly on the optical surface and 
focuses onto a focal “point” on the shroud where the breakdown of air starts. Ten laser rays are 
plotted for clarity, whereas two hundred rays are actually used in the computation. 

It is probably best to recount the phenomenon of optical breakdown, also known as the optical 
detonation, at this moment. This was put forward by Raizer"^’ in detail in his book. Essentially, 
a strong shock wave is generated in the region where laser radiation is absorbed by the plasma 
and heat is deposited very rapidly. This shock wave travels through the medium, heating and 
ionizing such that the medium becomes capable of absorbing more laser radiation. The laser 
energy is deposited in successive layers of a medium adjoining the front of the shock wave 
subjecting to the laser radiation. These layers themselves become energy sources which maintain 
the shock wave. Hence, the shock wave moves along the optical channel opposite to the laser 
beam, whereas during the initial stages it is maintain by the beam and does not decay. 

Figure 3 shows the temperature contours at an elapse time of 10 (xs where the optical 
breakdown is being fed by the laser energy and the plasma front grows. Notice the rays are 
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allowed to bend as the index of refraction varies with the expanding plasma. As the local 
electron density of a plasma front reaches the critical value such that plasma resonance occurs, 
the local plasma becomes opaque and the regular laser energy absorption moves to other area of 
the plasma front. This process is nonlinear and may explain the irregular shape of the plasma 
front. It is anticipated that laser beam totally reflects if the plasma resonance spreads to the entire 
plasma front. The phenomenon of plasma resonance disappears as the plasma front expands 
continuously, while the electron density decreases continuously. The computed optical 
breakdown phenomenon appears agreeing with that described by Raizer. Notice the spreading 
of the shock wave outside the optical channel results in gradual attenuation because there is no 
energy to maintain it. In the downbeam direction the shock wave is either weakly maintained by 
the laser radiation than in the upbeam direction, or is not maintained by the laser radiation during 
plasma resonance. 

Figure 3 shows the temperature contours at an elapse time of 20 p,s. At this time, the laser 
beam is turned off since the pulse width is fixed at 18ps, thus the shock wave is not being 
maintained by the laser radiation. The computational results show that a pulsed detonation wave 
usually takes place at approximately 10 to 20 ps after start command. The higher the pulse 
energy, the earlier the generation of a pulsed detonation wave. Part of this propagating 
detonation wave will hit the parabolic optic and reflect off the surface. The reflected and non- 
reflected portions would combine to make a propagating Mach stem. It should be noted that the 
phenomenon of a propagating Mach stem on the optical surface is very similar to that of a 
propagating Mach disc inside a bell nozzle during the start-up transient.''^’’''^^ This is because the 
aerospike optic geometry is an “inside-out” version to that of a contoured bell nozzle. 
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As the Mach stem moves past the end of shroud, the part that is attached to the inner-shroud 
starts to wrap around the outer-shroud. As the wrap around progresses, the shock wave weakens 
and grows like an ellipsoid. A forebody attaching “leg” would develop, as shown in the pressure 
contours of Fig. 5. Eventually, the “leg” leaves the nose and the ellipsoid expands to infinity. 
That completes the propulsion physics of the Laser Lightcraft for a single pulse of laser energy. 

Coupling coefficient is the generated thrust divided by the delivered average laser pulse 
energy. It is a performance measure particular to the pulsed laser propulsion at the air-breathing 
mode, and conceptually similar to the specific impulse (generated thrust divided by mass flow 
rate) of a chemical rocket. Figure 6 shows a comparison of the coupling coefficients between the 
model prediction and test measurement. The model predicted coupling coefficients agreed 
reasonably well with those of the measurement. At this writing, the effect of non-equilibrium 
radiation has not been added to the result, hence there is room for improvement once it is 
considered. 


Conclusion 

A computational plasma aerodynamics model has been developed to study the propulsion 
physics of an experimental Laser Lightcraft. The model development is based on a build block 
approach such that the model can be improved continuously with improved understanding of the 
physics. The model predicted laser propulsion physics such as the optical breakdown and 
detonation wave propagation agree well with those described in the literature. The model 
predicted coupling coefficient for a Model #200 Laser Lightcraft agreed reasonably well with 
those of the measurement. 
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Appendix 1. Air breakdown kinetics mechanism. 
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Fig. 1 The layout of a computational grid for the laser lightcraft. 



Fig. 2 Computed temperature contours and laser ray traces at 0.5 |is. 





Fig. 3 Computed temperature contours and laser ray traces at 10 jis. 





Fig. 4 Computed temperature contours at 20 ps. 





Fig. 5 Computed pressure contours at 105 fis. 
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Fig. 6 A comparison of the coupling coefficients. 




